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The behaviour of a "test" electromagnetic field in the background of an exact gravitational plane 
wave is investigated in the framework of Einstein's general relativity. We have expressed the general 
solution to the de Rham equations as a Fourier-like integral. In the general case we have reduced 
the problem to a set of ordinary differential equations and have explicitly written the solution in 
the case of linear polarization of the gravitational wave. We have expressed our results by means 
of Fermi Normal Coordinates (FNC), which define the proper reference frame of the laboratory. 
Moreover we have provided some gedanken experiments, showing that an external gravitational 
wave induces measurable effects of non tidal nature via electromagnetic interaction. Consequently 
it is not possible to eliminate gravitational effects on electromagnetic field, even in an arbitrarily 
small spatial region around an observer freely falling in the field of a gravitational wave. This is 
opposite to the case of mechanical interaction involving measurements of geodesic deviation effects. 
This behaviour is not in contrast with the principle of equivalence, which applies to arbitrarily small 
C^i . region of both space and time. 

<N ; 
X> : 

D ■ I. INTRODUCTION 

\Q ' The behaviour of a "test" electromagnetic field (i.e. an electromagnetic field whose stress-energy tensor does not 
affect the curvature of the underlying space-time) in a gravitational wave background has been widely studied in 
the framework of Einstein's general relativity. This topic is very important to conceive possible further experimental 
verification of general relativity and also to better understand the principles underlying the theory itself. A great 
deal of efforts have been aimed at solving Maxwell (or de Rham) equations to first order in the gravitational wave 
amplitude [p|-|l2| or, at most, to second order under geometrical optics limit |ll|. However, a general solution within 

^v^j . the framework of the full theory of general relativity would highlight the main features of a free electromagnetic 
field in radiative curved space-time. Indeed any possible ambiguity arising from approximation procedures could 



be circumvented 1 1 1 12[] . We point out that such ambiguities exist not only when strong gravitational waves are 
concerned, but also may appear for weak gravitational field s as well. An example of this is the question related 



q \ to photon creation [see discussion in Sec. VI after Eq. (6.27)] . Therefore any result obtained in the framework of 



H ' linearized general relativistic theory has to be handle with care due to the non linear character of the full Einstein 
theory. The linearized theory is not general relativity but minkowski flat space-time plus a small perturbation. 

To solve the problem one can take advantage of the freedom in the choice of a particular system of coordinates 
in which writing the equations. Obviously it is convenient to choose the reference frame in which the equations are 
easier to be solved. This is usually accomplished when the metric tensor is expressed in its simplest form. This system 
attached to the wave will be hereinafter referred to as Wave Reference Frame (WRF). In the limit of linearized gravity, 
it reduces to the usual Transverse Traceless (TT) gauge |Q. Recently a solution in the WRF has been obtained p5[ . 
However, such a system is not the one where measurements are performed (laboratory frame). Therefore, in order to 
obtain results that have a direct interpretation from the physics point of view, one should express the electromagnetic 
field variables in a reference frame attached to an observer. The most natural way to construct such a frame is to 
consider the observer freely falling in the field of the gravitational wave. Such coordinates are the well-known Fermi 
Normal Coordinates (FNC) 

The aim of this paper is to obtain the solution of this problem in FNC and employ this solution to predict a few 
new physical effects. To this purpose we have re-obtained the solution in WRF in a slightly different way, expressing 
it as a Fourier-like integral. This form is more suitable for applications to concrete physical situations. Moreover 
we have obtained the transformation rules between WRF and FNC which hold true in each point of the region of 
space-time where the plane propagation of the exact gravitational wave is valid. 

One of the most important results is that in FNC it is possible to infer the presence of a gravitational wave over an 
arbitrarily small spatial region in the neighbourhood of the origin. In other words there are both tidal and non-tidal 
effects in spatial coordinates. Therefore, near the origin, d'Alembertian operator applied to the electromagnetic field 
is non- vanishing, but rather proportional to the order of magnitude of the electromagnetic field. This is not surprising 
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because Maxwell second-order equations involve the Riemann tensor J14:|Jl 7| ] . Consequently the solution to Maxwell 
equations can not be written as a monochromatic plane wave and no photon is created by the interaction. We stress 
that this behaviour is not in contrast with the principle of equivalence, which applies to arbitrarily small region of 
space-time . This topic is more deeply investigated in jy! . 

The paper is organized as follows. In Sees. II and III we perform calculation in WRF. In Sec. IV we provide 
the transformation rules connecting WRF and FNC. In Sec. V we express the four-vector potential in the latter 
frame. Finally, in Sec. VI the theory is applied to a case that is theoretically interesting and opens up some gedanken 
experiments. 



II. SOLUTION TO THE FREE DE RHAM EQUATIONS IN WRF 



Let us call the coordinates in WRF. According to @,|9 20 an exact plane gravitational wave propagating along 
the y 3 -axis is described by the following line element 



ds 2 



V) 2 + (y 3 ) 2 + 7V-yW) 2 
My 3 - V°)(dy 2 ) 2 + 2S(y 3 - y^dy'dy 2 



(2.1) 



We do not make any assumption about the actual form or amplitude of the three quantities V, <2, and S but simply 
require that the metric tensor is an exact solution to the vacuum Einstein equations. Moreover, were the curvature 
small, the line element (2.1) would just describe a weak plane waves in a TT reference frame [Q. Since metric 
coefficients depend on the coordinates only through y 3 —y we are naturally led to perform the following transformation: 

Y°=y 3 -y° 



Y 1 
Y 2 
Y 3 



y 

y 3 + y° 



(2.2) 



Since the wave keeps plane in the new set of coordinates, the new reference frame Y^ will still be referred to as WRF. 
The covariant components of the metric tensor become: 



( 

V(Y°) 

S(Y°) 

h o 



o 





S(Y° 
Q(Y°) 
/ 



(2.3) 



In order to have a real space-time, the following inequality must hold true: 

1 



.9 



(TQ-S 2 ) < 



(2.4) 



First of all we consider the general form of de Rham equations under Lorentz condition, describing the propagation 
of free electromagnetic field in empty curved space-time (we stress again that the electromagnetic field strength is 
considered so small as to neglect its effects on the curvature of space-time): 



A a ' v , v = 
A».„ = 



(2.5) 
(2.6) 



By substituting the metric tensor ( [2.3] ) in (2 1 5)-(2 ; 6) we achieve (we use convention and notation as in |14[| except for 
p, q, r, and s which take values 1 and 2): 



g^A\. v + 2 g^g qp d 3 A p + g ap g pq d r A° - 2 5% g^g rq d p A r 
1 



■S^g^g^g^A 3 



2 dY° d3A 



0, 



(2.7) 



with Lorentz condition: 



dY° 



A° = 0: 



(2.8) 
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dot means derivation with respect to the argument. As is well known in classical electrodynamics (e.g. pi]] ), even for 
potential satisfying Lorentz condition there remains the possibility to perform a restricted gauge transformation such 
that the resulting new potential still satisfies the Lorentz condition. 

Our aim is to solve the system of equations (2.7) taking into account condition (2.8); in order to simplify them, 
we choose the particular restricted Lorentz gauge in which .4° = 0. We will show that this choice is always po ssible, 
without loss of general ity (analogously to the classical charge-free case |2l]]). In fact, if A'^ satisfies Eq. (2/7) with 
the Lorentz condition (2.8), then also A p + A ;A1 does, provided that: 



g \ v ,oi + Z ^tj; A i3 — u, 



(2.9) 



wher e A is a scalar function. Setting A = 0, then A 3 = \A' a . If this last condition holds, then, by substitution in 
Eq. (2/7) for a = 0, it follows that A indeed does satisfy: 



^W2^I^A 3 



= 0. 



(2.10) 



' ,3 

In conclusion, for the charge-free ca se, it is al way s possible to set A = 0. 

Because of this choice both Eqs. (2.7) and (2.8) turn into simpler ones. Namely: 

gTA"^ + 2g aq g qp d 3 A p - 2 6% g qp g rq d p A r + 2 d 3 A a = 



and: 



d k A k = 



(2.11) 



(2.12) 



The advantage of having replaced the old variables y a with Y a is that in Eqs. ( 2.11 ) there is no second-order 
derivative with respect to either Y° or Y 3 . This fact, together with the dependence of t he co efficients of the metric- 
tensor components only on Y , make it possible to turn the partial derivatives system (2.11) into a simple ordinary 
one. If we write the solution as a Fourier integral 



A a = / d 3 Ae[^ YJ )]a tt (y ,A 1 ,A 2 ,A 3 ) 



(2.13) 



and substitute it in the first two equations of system ( 2.11 ), we get rid of the dependence on Y 1 , Y 2 and Y 3 and 
obtain the following ordinary system for Fourier- coefficients a 1 and a 2 : 



4*A 3 ^ - (A : 2 + A 2 2 ) a s + f pq \ P \ q a s + 2zA 3 g^g pq aP + iA 3 d ^ g) a s = 

where we have set f qp = 8 qp — g qv . The above system of equations can still be simplified by setting 

a s = jB s 

with 
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(-<?)' 



exp 



{x 2 + \ 2 ) r" - f^ q x P \ q 

4iA 3 



(2.14) 



(2.15) 



(2.16) 



where 1 - )p9 is any primitive of f pq . Equations (2.14) yield the following system of equations for the unknown 
quantities B s : 



dY° 2 p 



where we have set: 

M\. = g sq g qp 

before discussing the solutions of system ( [2.17| ) we consider last equation of the system (2.11): 



(2.17) 



(2.18) 
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4*A 3 ^ - {X 2 + X 2 ) a 3 + f pq X p X q a 3 - 2iX p g™g qs a s + tX 3 g) a 3 = (2.19) 

The solution of this equation comes directly from Lorentz condition: 

a° = 0; a a = 1 B 3 ; B 3 = -^- (2.20) 

Xr 



Now only Eq. ([2.17D needs to be solved: once B 1 and B 2 have been obtained, it is trivial to achieve A a by means 



of Eqs. J2.15D , (p.20D and ( E.13| ). The solution to Eq. ( |2.17|) can be more or less difficult to be deter mined depending 
on the actual form of V(Y {> ), Q(Y°), and S(Y°). In the general case, an analytical solution to Eqs. ( 2.17 ) could not 



be available. A major advantage of this system is its n on-dependence on therefore a numerical solution can be 



easily implemented. The general solution to Eqs. ( 2.17 ) can be inferred by multiplying the numerical solution by an 
arbitrary function of A^. 

The results of this Section are in agreement with those obtained in [[L5|| . Yet, we have used a different formalism 
which is more suitable for ap plication to concrete physical situations. For instance, by choosing properly the Fourier- 



like coefficients of Eq. ( 2.13 ), any possible particular boundary or initial condition can be obtained. 

In next Section we will show that it is possible to solve exactly such a system for the case of a linearly polarized 
gravitational wave. 

III. LINEAR POLARIZATION 

An exact plane gravitational wave is linearly polarized when the metric tensor can be put into diagonal form through 
a fixed rotation in the Y x -Y 2 plane. 

Consequently, without losing generality, we choose the Y^ reference frame in such a way that the only non-vanishing 
components of the metric tensor are: 

_ 1 

.903 — .930 — 2 

gu=V{Y°), 522 = Q(y°). (3.1) 

If we set V(Y°) = L 2 (Y°) exp(2/3(r )) and Q(Y°) = L 2 (Y°) exp(-2/3(r )), L(Y°) and p(Y°) being arbitrary 
functions, then the condition for the gravitational field equations to be satisfied becomes Q: 



d 2 L ( df3 y 
d(Y ) 2 + \dY° ) L ~ 



L = (3.2) 



One can easily see that, were f3 small, the metric would reduce to that of a usual weak gravitational wave in linearized 
theory |Q. According to L(Y°) and f3(Y°) are referred to as background and w ave factors. 

We now proceed to calculate the exact solution to the system of equations ( [2.17j ). By using the background and 
wave factors instead of V and Q, it reads: 

B 1 (Y°) = (3.3) 
B 2 (Y°) = 0. (3.4) 



dB 1 


1 


( dL 


f L(Y°) 


d(3 


dY° 


V L 


\dY° 


dY° 


dB 2 


1 


( dL 


-L(Y°) 


dp 


dY° 


V L 


(dY° 


dY° 



If we set: 



V + = P+\ lg(L 2 ), V- =P- l -\g{L 2 ) (3.5) 



after some manipulation, we obtain: 



dB 1 
B 1 
dB 2 



-dV+ (3.6) 
dV-. (3.7) 
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These equations can be easily integrated obtaining: 

B 1 (Y°)=b 1 e- v+ (3.8) 
B 2 (Y°)=b 2 e v ~ (3.9) 

where b 1 and b 2 are constants. Substitution of Eqs. fl3.8|) and fl3.9|) in Eqs. (2.20) yields: 



B A = --^b'e- v --^b'e v . (3.10) 
A3 A3 



As for Eq. ( 2.16 ), we obtain 



where V + and V are functions defined as follows: 

dV ± (w) 



dw 



(3.12) 



IV. FNC FOR AN EXACT LINEARLY POLARIZED WAVE 

The results we have obtained in last section are valid in the WRF; though they do not correspond to any measurable 
quantity. In fact, the proper reference frame where an observer may execute a measurement is the FNC |l^ , |l6| . 
Therefore any tensor expressed in WRF must be written in FNC. To this aim, this section is devoted to linking the 
WRF to the laboratory reference frame (FNC); this is accomplished by following the procedure outlined in |lj,|l6|. 
First, we consider an observer moving along a time-like geodesic in WRF; let q(r) be its world line, r its proper time, 
and /(o) its four-velocity. Thus we can write: 

/(o)- dT ■ (4.1) 
To achieve the actual form of the geodesic the following equations need to be solved: 



dm 



(0) 



W<V(0) = "I (4-3) 

The only non-vanishing Christoffel symbols are: 

, dV+ 2 dV- 

1 01 — «" " 1 n 



dY° vz 3Y° 



F 3 d _ ( p 2V+(Y°)\ p3 ^_ ( -2V-(Y°)\ (44) 

1 11 - ay° \ ) 1 22 - oy° \ e ) 1 ' 

After straightforward calculation we obtain: 

/(°o)(t)=^o) (4-5) 

/ ( 1 0) (r)=^ 0) e- 2V+ Ko,^ (4.6) 

/( 2 o)(r) = ^o)e 2V_( ^ T) (4.7) 

fh(r) = -^e-^Ko^) .W^r)^ (4 . 8) 

V(0) ^(0) 

where </3" ^ are four constants. The geodesic parameterization may be achieved by performing the following integration: 

q<*( T )=q a (0) + f f^(r')dr'. (4.9) 
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Once /(o) has been determined, which is the time-like vector of the orthonormal tetrad carried by the observer, we 
determine the other three orthonormal space-like vectors /(j), to complete the observer's orthonormal tetrad |]l^ , |l6| . 
These vectors have to be parallel transported along q(r). Therefore the following set of equations must hold true: 



^ + r^/ ( yf 0) = o. 



The solution can be written as: 



1 } , 

tfi)rfo) p -2V+(y° OT r) 
-,0 ,„2 



/2 ( T) = WW eaVftpfo r) + 2 e V-(^ 0) r) 



V(0) 

,3 r v _ y W ( ^(0) )2 _-2V+(^ 0) r) _ y(i)(y(0)) 2 2V-( V ° (0) r) 

hi){ ) Ko)) 2 ~ Ko)) 2 

_ 2 y>y ) e - V+ « 0) r) _ 2 yWo) e V-(^ r) 3 

where y?U are twelve constants. Besides, orthonormality conditions imply: 

The conditions given in Eqs. (4.5) and ( |4 . 1 5[ ) hold true provided that: 

4) P(o) + V?o V( 0) = o, 

4) + 4) 4) + 2 Ptf) + 2 4) 4) = 2 



(4.10) 

(4.11) 
(4.12) 

(4.13) 



(4.14) 



(4.15) 



(4.16) 



Next step is construction of a space-like geodesic, originating from an arbitrary point on the observer world line. 
Let a be the space-like geodesic parameter; therefore Y(a, r) is the point parameterized by the cr-variable on the 
space-like geodesic whose emanation point is q(r). Let f a — be the tangent vector in Y(a,r) to the space-like 



geodesic; the geodesic equation reads: 

da ^ J J 

We seek out the solution to the above equation satisfying the following conditions at a = 0: 



df c 



F(a = 0)=a 3 f^(r) 
F> = 0) = g"(r), 



(4.17) 



(4.18) 
(4.19) 



where a 3 are the direction cosines of the geodesic at q (r). In order to achieve / we perform a calculation, which is 
formally very similar to the one we made to solve Eq. (4.2); then we obtain Y a (a, r) by direct integration: 

Y a (a,T) = q a (T) + T f a da' (4.20) 
Jo 

A natural way to identify an event in the observer's reference frame is taking the four numbers [fl4|| : 

(4.21) 



x° = r 



Substitution of these identities in Eq. (4.2C) yields the desired coordinate transformation: 



G 



yO 
Y 1 



Y' = 



^U^) -V+(^ ,x°) 

( i P(k) x ) i- 



V+(^ 0) x° + ^ k) x k )-V + (^) 



V + (^ a) x° + ^ k) x k )-V + (0) +q 1 (0) 



¥>(0) 



V-(v i0) x + ^ k) x k )-V-(v° i0) x ) 



V-(v io) x + i p {k) x k )-V-(0) +q 2 (0) 



Y A = - 



H)) 2 

Ko)) 2 



V + (^ 0) x°+^ k) x k )-V+(0) 
i „k ,„i 



(^(fc)^^ 2 c 2V+(y° m x°) , ■> ""'V" ^''n ,V (/,'„■", 

(y>° x fe ) 2 



(y(fc)* fc ) 2 



^(fe)^ n 0) 

-2V-(y° 0) x°) g ^(fc) xfc ^(0) _-V-(^ ,»°) 
' h 

^(fc)* V(0) 



y-(^ 0) o: + ^ fe) x fe )~y-(o) 

V-(^ )*°+*'(*)^)-^-(^o)^ ) 



(V( fc) ^ fe ) 2 (^( fe) ^) 2 



In the framework of weak-gravitational-field approximation a similar result was found in ]l6[] . 
The metric tensor components in FNC could be written through the usual relation 



dY a dY? 



dx^ dx 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



Yet, for our purposes, it is sufficient to know that 7 M „ = r/^ + O (|a;| 2 ) and we do not need to calculate them explicitly. 



It can also be checked that r^o = O (\x\ 



Although coordinate transformation ( 4.22 )— ( 4.25 ) was determined in the case of a diagonal metric te nsor in WRF, 
the expression we obtained holds true for a generic linearly polarized plane wave. Consequently Eq. (3.1) does not 
lead to any loss in generality. 



V. FOUR-VECTOR POTENTIAL IN FNC 



The knowledge of the transformation rules given by Eqs. (4.22)— (4.25) allows one to obtain — via direct differenti- 
ation — the expressions of needed to calculate covariant components of a generic tensor in the FN C. T he refor e 
t he co v arian t c omp o nents of the four-vector potential in the WRF are to be determined. From Eqs. ( |2.3| ), ( [2.13 ), 
(|t|) , (|2~20| ), Q-(|T|) we get: 



ao = —- 



I- (\ 1 b 1 e~ v+{Y ^+X 2 b 2 e v ^ Y "A e -^ + O-Vv-(Y 
2 V / 



V2A 3 

a 1 =V2b 1 e^ v+ ^ +v ~ ( - Y °^e 



(>irv + (y")+(>2)"v~(y") 

)) g 4iA 3 



- tv o\\ (Ai) z v+(y")+(A 2 )^y-(y") 

4iA, 



« 2 

a 3 = 



4»A 3 



(5.1) 

(5.2) 

(5.3) 
(5.4) 



By setting as the covariant components of the four-vector potential in the FNC, they can be achieved by: 



" ~ dx» A ° 



(5.5) 



The components are written as: 
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where 



a 



b l e -i(3V+(0)-V-(0)) 



„2 



CO) ^ 



V2A 3 

+ V2^ ) [l> + (<A))e V+(W (^ + (0)-y+(0 o ))-e- v+ ^) 



fe 2 e -I(V+(0)-3V-(0)) 



,V-(0o) 



b 1 e K v+ w+v- w) 



6 2 e -i(V+(0)+V-(0)) 



(5.6) 



(5.7) 



and 



.1 J* 



~k 

„2 „,fc 



pg. ) 6 1 e-*< 8V+ W>- v ~W» 



y V2A3 



& 2 e -i(V+(^)-3V-(^)) 



(^° fe) a; fc ) 2 



* = ^ ( A: - 2A 3 §^ ] D+ + ^ {k) x k e~ v -^) [ x 2 - 2A, — • | D 



A- U 0) V+ + <p 2 {0) V-) + A 3 (C+ + C~) + A 3 ^^ fc + Aa^o)* + \ 3 q 3 (0) 
¥>(0) v ' 

(x 1 fv+(<f>) + (x 2 ) 2 v-(<i>) 



In the previous equations we set 



4A, 



00 = f° )X°. 



As for Eq. (5.9), the following functions were introduced: 



<fi° {k) x k 
(^( fc )^ fc ) 2 

h 



!_ e 2V+(0) D + 



V ± (<P)-V ± (ct> ) 



f° {k) x k 



<P(0) 



v + = L f l [v+(4>)-v+(p)], 

V(0) 
2 

p - = ^ [v-(<f>)-v-(o)]. 



(5.8) 



(5.9) 
(5.10) 

(5.11) 
(5.12) 
(5.13) 
(5.14) 
(5.15) 
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VI. APPLICATIONS AND DISCUSSION 



In order to emphasize geometrical effects over cinematic ones we are naturally led to choose the coefficients as 
follows: 



u 

4) = Di i 



(6.1) 



where 



Di k D i k = si r 



(6.2) 

With this choice the observer is assumed to be at rest in the reference frame with given orientation with resp ect to 
in q(r). Euler angles of rotation matrix D l k determine the orientation. One can easily see that the conditions (4.1£) 
are met. 

In order to obtain the relationship between Fourier coefficients A and the usual flat space-time wave-vector, we set 
L = 1 and /? = in formula ([5~9|): 



(Ar) 2 + (A 2 ) 2 n3 
4^ Dk 



'" + A 3 + 



(A!) 2 + (A 2 ) 2 \ „ 



4 A, 



In a flat space-time, is usually defined as 

* = * ± = kjX j ±kx° k = ^kf + kg + kg 

From the above equations, by taking into account the orthogonality of matrix £)* ■ we obtain: 

A, 



As = A 



± D 3 J kj ± k 



The integral of Eq. (5.6) is expressed in terms of k as a sum of two terms: 



-4, 



d 6 k 



e o„ + e a. 



where a* = a M [A ± (fc), x Q ; &±(fc)] with a M given by Eqs. (5.7) and (5.8). 
In order to better understand the behaviour of the solution, let us suppose: 



In this case 



while 



k^±k 



= k 1 x 1 e v+ ^D+ + k 2 x 2 e- v ^-^D- + (x 3 + x°) 

{ktfV+jx 3 -x°) + (k 2 ) 2 V-(x 3 -x°) k 3 ±k (r+ , r . 

2(k 3 ±k) + 2 [L +L 



(6.3) 



(6.4) 



(6.5) 
(6.6) 

(6.7) 
(6.8) 



(6.9) 



where 
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± _ V ± (x 3 -x°) -V ± (-x°) 
x 3 

V ± = 



C+ 



c~ = 



(x^f 



{x 2 ) 2 



1 _ e 2V+(x 3 -x°) D + 



l- e -2V-(x 3 -x°) D - 



(6.10) 



As an example we choose the constants 6± so that there is a static magnetic field along the x 3 direction in absence 
of gravitational wave. This is mathematically accomplished by setting b r ± = b± <5[, where 



b+ = b_ = 



iB 



5{k 3 -k')6(k 1 )6(k 2 ). 



2V2k 

By performing the integration and taking the limit k' — > 0, we achieve: 



(6.11) 



Ao = _^! BojD - e -V-(-x°) ( e V+(-xO) e -i 



[3V+(x 3 -x°)-V-(x 3 -x )] 



V + (-x u )e 



A 2 = 



0^V+(-x°) / T/ +^3 



-v+(-x°) 



(y+(x 3 -x )-V + (-x ))-e- 
i[v+(x 3 -x°)+v + (-x )] e i[v-(x 3 -x°)-v-(-x°)] 



f; 2 



[V+(x 3 -x a )+V-(x 3 -x 



°)]J 



1 3 = ^f!^ 0j D- e V + (-^) e -V-(-x°) j e -i[3V+(x 3 -x°)-V-(x 3 -x )] _ jD+e i[V+(x 3 -x )+V-(x 3 -x°)]| _ 



(6.12) 

(6.13) 
(6.14) 



The above expressions hold true everywhere. In order to show their good behaviour in the neighborhood of the origin 
we perform a power expansion in x 3 up to \x\ 2 terms. We obtain: 



A =x 1 x 2 B V + {-x°)e-i 

.r 3 

1 



A 1 = x z B 
A 2 = 



i [V+(-x°)+V-(-x )] 
X " ' 0\ i 3 3 • . \ 



2 V+C-^ + ^a: 3 



-[V+(-i°)-2V-(-i )] 



A 3 = -x 1 x 2 B V+(-x°)e-^ v+ (-^- 3V ^-^ 
The electromagnetic tensor field components can be obtained from the usual relation: 

F —A 



_ A 



(6.15) 

(6.16) 

(6.17) 
(6.18) 

(6.19) 



Although it would be possible to calculate F^ v everywhere, we shall limit ourselves near the origin. Therefore, up to 
linear terms in x k , we achieve: 



F m = x 2 Bo 



V + (-x°)-2V-(-x°) 
x 2 B Q V+(-x°)e-t 



e 

Mv+(-*°)+v-(-x )] 



[V+(-x°)-2V-(-x )] 



F 02 = -x 1 Bo V+(-x°) e~i [v + (-°)+v-(-x°)] 

^03 =0 

Fu = —Bo 

F 23 - -x 1 Bo V+(-x°) e-i [v + (-°)-3 v-(-x°)] 



l-^V+(-x a ) + ^x 3 V-(-x°) 



-[V+(-x°)-2 V-(-x )] 



F 31 =x 2 B 



\v + {-x°) + \v-{-x«) 



-[V+(-x°)-2V-(-x )] 



+ x 2 Bo V+(-x°) e~i [v+(-x°)-3V-(-x°)] 



(6.20) 
(6.21) 
(6.22) 

(6.23) 
(6.24) 

(6.25) 



We notice that in the laboratory frame the effect of a gravitational wave on an electromagnetic field is not only of 
tidal nature. In fact when x k — > the the only non-vanishing component becomes: 
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Fv. 



-B{x° 



B e 



-{V+{-x°)-2V- (~x a )\ 



(6.26) 



Two important conclusions may be drawn from this rath er sim p le cas e. First, the solution can not be interpreted as a 
photon creation. In fact, the solution described by Eqs. ( |3.2C )-( |6.25|) has non-vanishing d'Alembertian, proportional 
to the field itself. Moreover direct computation of the divergence of Pointing vector at the origin reads: 



rriOk 

1 .k — 



dx° I 8 



Bp e -2[V + (-x°)-2V- (-2 



(6.27) 



It is the time derivative of a function assuming finite values. Its time average vanishes: in fact, if the gravitational 
wave is a periodic function, the time average of Eq. (6.27) over the period T is zero; for a generic non-periodic wave, 
the time average vanishes over an interval of time that is longer than the signal duration (see for instance JljJ, §34). 
Therefore there is no net flux of electromagnetic radiation across any close surface. This implies the absence of any 
radiation field. This conclusion is not surprising since the solution still keeps the tensor form of a static magnetic 
field. 

This is a good example in which the calculation in the framework of the full theory of general relativity settles 
down an open question of the linearized theory. In fact, by using the linear approximation, the time average of 
the divergence of the the Pointing vector at the origin vanishes (taking only first order terms). The first non-zero 
contribution is due to second order terms, which, however, are to be neglected in a linearized theory. One could be 
led to think that quadratic terms could possibly give rise to a non- vanishing time average of T k . In the full theory, 
Eq. (6.27) immediately shows that this is not the case. 

Another interesting feature is the possibility for an observer in FNC to prove — at least in principle — the presence 
of a gravitational wave, by performing an experiment involving electromagnetic interaction. The non-tidal nature 
of the interaction allows one to asses the presence of a gravitational wave irrespective of the apparatus size. This is 
quite different to mechanical detectors, involving geodesic deviation measurements. In the following we will propose 
a few gedanken experiments aimed at showing this peculiarity. 



A. Motion of a particle on the x = plane 

As a first example we consider a charged particle moving on the a: 3 = plane. The exact equations of motion are 
given by pi: 



(6.28) 



As = 0(x k ), we can assume that near the origin the leading term is the non-tidal one [see Eq. ( |6.26 )]. Therefore 
equations of motion become (7= (1 — f3 2 )~i , (3 = -, v is the speed of the particle): 



du° 
ds 
du 3 
ds 
du 1 
dx° 
du^ 
dx° 



u 3 = 



x° = 7 s, 
x 3 = 0, 



-~B(x°)u 2 , 
\B{x°)u\ 



(6.29) 
(6.30) 
(6.31) 
(6.32) 



where £ = 7 m c 2 is the energy of the particle. We see that the speed of the particle is a constant of motion. Setting 
B dx° = dw we get: 



u 1 = 7/3 ci cos ip(x°) + 7/3 C2 smip(x°), 
u 2 = 7/3 ci smip(x°) — 7/3 C2 cos tp(x ), 



B{Z)d£ + w a 



Ci 2 + Co 2 = 1 



(6.33) 
(6.34) 

(6.35) 



the particle performs a non closed orbit around the origin, with a variable period T(x°) given implicitly by: 
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r x°+cT(x°) 

2vr= / <p{€) d i (6-36) 



By measurement of the time variation of the period (which is constant in flat space-time) one can infer the presence 
of a gravitational wave. 



B. Induced e.m.f. in conducting circuits 



As a second example we consider a conducting ring with resistance R, lying in the x 3 = plane with its center at 
the origin (in this application SI units are used). Assuming the ring sufficiently small, we can neglect all the terms in 
the electromagnetic tensor except the one given by Eq. (6.26). Therefore Faraday law causes a current to flow. One 
has: 



/(f) 



5 dB{t) 
R dt 



(6.37) 



where S is the ring surface and x° = ct. 

In general one expects the same kind of effect in RLC circuits. Let us consider a series circuit with a resistance R, 
an inductance C, and a capacitance C under the same assumption as before. In this case the equation of motion for 
the charge on the condenser plates is: 



Q + —Q + uj 2 Q = -^B 



(6.38) 



where To = uj q = -J^, and in this case the dot means derivation with respect to t. We notice that in flat space- 
time, where B — 0, there is no current flow ing. The effect of the gravitational wave is a current flowing through the 
circuit. In general the solution to Eq. ( |6.38 ) is given by: 



S(t) 



iuj B{oj) 



i — 

TO 



(6.39) 



(f(t) 



i 

/2rr 



In order to study the main characteristics of the current we assume the gravitational wave to have a typical frequency 
g . In this way the Fourier components of the magnetic field are approximately given by: 



B(u) 



2irB \ 5{w) - -f [5(lu - tog) - 5(w + u g )] 



(6.40) 



Substitution in Eq. ( 5.3S ) gives for the flowing current: 

2 



lit) 



S Bo ^ g fo 



2SB Q 



fo Q cos u t 



(6.41) 



where Q = 



is the quality factor of the circuit. Therefore it is possible, at least in principle, to detect a 



gravitational wave by measuring a current. It is noticeable that the effect is not affected by the size of the device, but 
only by the circuit parameters as it was already known in the framework of linear approximation (e.g. |p2| , p3| ). 



C. Paschen Back effect 



It is well known that the spectrum of Alkali atoms embedded in a static and homogeneous magnetic field shows 
frequency shifts. We are here interested in magnetic fields strong enough to let us neglect spin-orbit terms. Besides 
we assume that the system is so close to the origin that we do not have to take into account tidal terms (i.e. we can use 
the usual quantum mechanics in flat space-time). Therefore under these assumptions we start from the Hamiltonian 
describing the interaction with an external electromagnetic field in the non relativistic regime (e.g. p4|). Substituting 
the value of the four-vector potential given in Eqs. ( 3.15| )-( |6.18 ), and neglecting tidal effects we get: 
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H = n + sn (6.42) 

no = ^ + V(r) (6.43) 

SH = --^- B(t) (L z + 2 S z ) (6.44) 
2 mc V / 

where e is the electron charge, m the electron mass, and V(r) the atomic potential energy. One could solve the 
problem using the standard time-dependent perturbation theory. However for typical frequencies of the gravitational 
wave much smaller than orbital ones (condition that is expected to hold true in nearly all circumstances), the effect 
results in a time variation of the frequency shifts of the flat space-time spectrum. Thus we get: 

(°) fJ>BB(t) \e\K 

Under this assumption gravitational wave causes the distance between two lines of the usual Paschen Back spectrum 
to change with time. 



D. Spin precession 



Let us consider a free particle of charge q and spin ~. Neglecting all its degrees of freedom except for the spin, the 
interaction with an external magnetic field can be described introducing the following Hamiltonian: 



ii = — — Bit) ■ S = — — Bit) S z 
mc mc 



(6.46) 



Once again we neglect tidal effects. This is true provided the particle is sufficiently near the origin. Therefore the 
magnetic field is given in Eq. (6.26). We obtain the evolution operator by solving Schrodiger equation. One has: 



at 



U 



-i I 0(t) Sz 



where 



eB 



[t + F(t)]; F(t)= f dt' 
Jo 



[ dt' 


\B(t) 1 
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L B 



2mc 

Let us suppose to prepare the particle in the following state (at time t = 0): 

la >= ci 1+ > +c 2 I- > 



(6.47) 



(6.48) 



(6.49) 



where |± > are the eigenkets of the spin operator along x 3 axis. Thus the state of the particle at time t is given by: 



\at >= U(t) \a >= ci e 



Let 



x+ >= 



1 

71 



> 



> +c 2 e 



>) 



Z0(t) 



> 



(6.50) 
(6.51) 



be the eigenket of the spin operator along x 1 axis with positive eigenvalue ]24| . As it is well known, the probability 
to find the particle in such a state at time t is given by 



p(t) = |< x + \at >|" = - [1 + 2 |d] |c 2 | cos(argc 2 -argci +2<f>(t))] (6.52) 
If the magnetic field were static (i.e. in flat space-time) this probability would periodically vanish when (n is integer): 

(6.53) 



mc 



t = TP = ("gci-argoOmc 

eB eB 



1 



2 MM 



2mc im 



eB n 



However, because of the presence of a gravitational wave this probability does not vanish in any longer, but takes 
values depending on n. Namely: 



V 



(4°)) 



1 — cos 



eB Q 



F(T^)) - V4MW-1 ( ^ F(Tf) 
/ \mc 



(6.54) 
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VII. CONCLUSIONS 



In this paper we have presented some interesting features of the behaviour of electromagnetic field in exact gravi- 
tational wave background by considering de Rham equations with no approximation. We have expressed our results 
in the laboratory frame (FNC) where any measurable quantity should be referred to. 

We have investigated a particular case in order to better understand the main features of the solution. In particular 
we have shown the appearance of non-tidal effects. Furthermore we have conceived some experimental setups which- 
in principle-could outline the differences between the response of an electromagnetic device and a mechanical one. In 
fact the latter is related to the geodesic deviation, while the former includes spatial point-like effects as well. This 
behaviour is not in contrast with the principle of equivalence, which applies to arbitrarily small region of space-time. 

It is also explicitly shown that in FNC the interaction does not create any photons. This is due to the fact that 
the d'Alembertian of the electromagnetic field does not vanish, being proportional to the field itself. 
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